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We present results for the confinement-deconfinement interface tension a^d of quenched QCD. They were 
obtained by applying Binder's histogram method to lattices of size L^ xL^x Lt for Lt = 2 and L = 8, 10, 12 and 14 
with Lz = 30 for L = 8 and L^ — 3L otherwise. The use of a multicanonical algorithm and cylindrical geometries 
have turned out to be crucial for the numerical studies. 
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1. Introduction 

At high temperature a phase transition oc- 
curs in QCD. In the quenched approximation (i.e. 
without any light quarks) this transition is of first 
order and separates a low temperature confined 
phase from a high temperature deconfined phase. 
The dynamics of a system which crosses the tran- 
sition temperature T^ (as e.g. in the early uni- 
verse or in heavy ion collisions) depends on the 
free energy 



Fed = ctcd A 



(1) 



of an interface of area A between regions of 
confined and deconfined matter. The interface 
tension acd = crcd Tc was investigated before 
in Monte Carlo simulations of lattice systems 
with Lt ~ 2 using various approaches (see |l|- 
Q). Lately these results have been questioned 
based on an application of Binder's histogram 
method to cubic spatial volumes L^ with L — 
6,8,10, and 12 (see ^j^)- However, these re- 
sults might have been plagued by interfacial in- 
teractions. Therefore, we present results using 
the same method but on asymmetric volumes 



(iz > L^ 
actions. 



Ly) thereby reducing these inter- 
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2. The Interfacial Free Energy 

We consider SU(3) pure gauge theory with the 
Wilson action S* on a cylindrical lattice of size 
LxXLyXLzX Lt with L^ — Ly — L and L^ > 2>L 
at the critical coupling Pc for Lt — 2. We use pe- 
riodic boundary conditions in the time direction 
and C— periodic boundary conditions Q in the 
spatial directions, i.e. 



Uf,{x + Lie,,t) = 


= U*{x,t), for 1 = 


^x,y,z (2) 


U^{x,t + Lt) -- 


= Uf,ix,t). 


(3) 



Because of the C— periodic boundary condi- 
tions the value of the Polyakov line i^L^x) = 

tr ( Y\t=i Uo(x, t) j will satisfy 

VtL{x + LiBi) = Vl\{x) ioY i^x,y,z (4) 

Therefore, no bulk configurations in either of the 
two deconfined phases that have nonvanishing 

imaginary part of Q.l = l/(-^^iz) X^s ^i(^) '^^'^ 
exist and the probability distribution PL{p)dp of 
p = Re VLl takes the form sketched in Fig.|l|. The 
system is most likely in either the one remain- 
ing deconfined phase corresponding to p'^^ or the 
confined phase at p^^>. When p is increased from 



n(l) 



bubbles of deconfined phase form. These 



configurations are suppressed by the interfacial 
free energy GcdA, where A is the surface area of 
the bubble. It grows until finally its surface is 
larger than the surface L^ of two planar inter- 
faces which devide the lattice into three parts 



Figure 1. Schematic probability distribution for 
the order parameter. The dotted line indicates 
the multicanonical distribution. 



Figure 2. Typical cuts through the lattice in 
the y — z— plane at values p'-^^ (first picture) and 
p{rain) ^ggcgud picturc) of the order parameter. 
The two phases are represented by white resp. 
shaded areas. 



as depicted in the second part of Fig. g. Since 
the interface area of the two planar interfaces is 
independent of p the probability P/, is constant 
in the region where their contributions dominate, 
i.e. around p('"™). Because of the C— periodic 
boundary conditions these interfaces always sep- 
arate a region in the confined phase from one in 
the deconfined phase that has Imi^h) = 0. Thus 
the corresponding configurations will be exponen- 
tially suppressed by the interfacial free energy of 
two confined-deconfined interfaces. Taking into 
account the capillary wave fluctuations of the in- 
terfaces as well as their translational degrees of 
freedom leads to additional power law corrections 



givmg 



PI 



ex L\ ■ L'^-^ ■ exp (-2CTcdi ) 



(5) 



for d— dimensional spatial volumes. This relation 
will be used to determine Gcd from the distribu- 
tions obtained on finite lattices. 

In order to calculate the probability distribu- 
tion Pl{p), one has to simulate the SU(3) pure 
gauge theory at the deconfinement phase transi- 
tion. But because of eq. a any standard local up- 
dating algorithm will have autocorrelation times 
tl which increase exponentially with L^ (" super- 
critical slowing down" ) . The use of the multi- 
canonical algorithm reduces this effect consider- 



ably. 

3. The Multicanonical Algorithm 

In order to overcome the supercritical slowing 
down, the multicanonical algorithm |l^Jll|l does 
not sample the configurations with the canonical 
Boltzmann weight 



\S) ex exp(/5S'), 



(6) 



where S — 1/3 ^p ir[/n is the Wilson action 

in four dimensions, 

weight 



but rather with a modified 



Vr{S)o,eMl3L{S)S 



«l(^)). 



(7) 



The coefficients a^ and (3l are chosen such that 
the probability Pl (not to be confused with the 
Boltzmann weights) is increased for all values of 
the action corresponding to the region between 
p(^^ and p^^\ as shown schematically in Fig. |l|. 
Details are described in Q where the efflciency 
of the multicanonical algorithm for SU(3) pure 
gauge theory has been demonstrated. We apply 
the algorithm to the determination of the inter- 
facial free energy. 

4. Numerical Results 

We have determined the probability distribu- 
tions for Lt = 2 and the spatial volumes L^ x L^ 



Figures. Eef^LCz) for a 14^x42x2 lattice. The 
dotted lines indicate the bulk expectation values 
of the two phases. 



Figure 4. Measured probability distribution for 
the order parameter. 



with L = 8, 10, 12, and 14 and L,_ = 30 for i == 8 
and Lz ^ 3L otherwise. Fig. g shows the real 
part oiflLiz) ee l/L^ J2x,y^L{x,y, z) for a typi- 
cal configuration close to p(™™) on a 14^ x 42 x 2 
lattice. As expected from section |[ one can iden- 
tify two interfaces between the confined phase and 
the deconfined phase. The imaginary part of fii, 
is always zero. In Fig. ^ the resulting probability 
distributions are shown. In contrast to the dis- 
tributions for cubic volumes {L values as before) 
they all have a region of constant probability in 
between the two peaks. This supports the sce- 
nario developped in section ||. 

In order to extract the interface tension we 
evaluate the quantities 



of 1/L'^. Their intercept with the y-axis is CTc 



^(1) 



(2) 



We extract Fj^ ' and F\^ ' from the probability 
distributions of Fig. Q. The results are plotted in 
Fig. |5| together with the corresponding linear fits. 
For the interface tension we get the value 



^=0.10(1) 



(10) 



It agrees within errors with the value obtained by 
replacing the Polyakov line by the action density. 
The agreement with [||-|[ is good while [Q quotes 
a slightly higher value. Still the discrepency be- 
tween these results and [pUq] which used Binder's 
histogram method for cubic volumes is reduced 
considerably and can thus be attributed mainly 
to interfacial interactions. 
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where Pl 



m 



X _ 1 
——(max) 

one expects Fl oc 



pr 



(8) 



(9) 



Note that 



\/L^I?. According to 
eq. H both quantities should be linear functions 
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